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In the first part of the talk, we briefly review the problem of parton-energy 
loss and thcrmalization at the Relativistic Heavy-Ion Collider and discuss how 
heavy quarks (charm and bottom) can help to resolve the existing experimental 
and theoretical puzzles. The second part of the talk is devoted to the properties 
of heavy quarkonia in the (strongly interacting) Quark-Gluon Plasma (sQGP) 
and their consequences for observables in heavy-ion collisions. 
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1. Introduction 

Experiments at the Relativistic Heavy-Ion Collider (RHIC) suggest that the 
matter created in semi-/central Au-Au collisions at v /SAw=200GeV r consti- 
tutes an equilibrated, strongly interacting medium at high energy density, 
well above the critical one at the expected phase boundary, e c ~ 1 GeV/fm 3 
(cf. the experimental assessment papers 1 ' 2 and references therein). This 
conclusion is essentially based on three evidences: (i) at low transverse 
momenta, p<r<2-3GcV (comprising ^99% of the produced hadrons), ideal 
relativistic hydrodynamics describes well the single-particle spectra and 
their azimuthal asymmetry, V2(pt)] the underlying collective expansion 
follows from kinetic equilibration times of around ro~0.5fm/c, implying 
rapid (local) thermalization of the system at initial energy densities of 
about eo— 30 GeV/fm 3 . (ii) at high momenta, pr>5-6GeV, the production 
of hadrons is suppressed by a factor of 4-5 relative to elementary p-p (or 
d-Au) collisions. This has been interpreted as energy loss of energetic par- 
tons traveling through an almost opaque high gluon-density medium, (iii) 
at intermediate momenta, pt— 2-6GeV, surprisingly large baryon-to-meson 
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ratios, as well as an empirical scaling of the hadron V2 according to its 
constituent quark content, has been interpreted as quark coalescence as the 
prevalent hadronization mechanism in this regime, and thus as partonic 
degrees of freedom being active. The microscopic mechanisms underlying 
thermalization and energy loss, however, are not well understood yet. In 
fact, the jet quenching of light partons is so strong that high-p<r particle 
emission appears to be mostly limited to the surface of the fireball, ren- 
dering a precise determination of the corresponding transport coefficient 
(and maximal energy density) difficult. While earlier approaches assumed 
the dominance of radiative energy loss via medium-induced gluon radia- 
tion (which prevails in the limit of high jet energies) cf. 3 , it has been 
realized subsequently that elastic (2<->2) scattering processes cannot be ne- 
glected at currently accessible pt<10-15GgV at RHIC. Moreover, none 
of the evidences (i)-(iii) is directly connected to the fundamental proper- 
tics distinguishing the Quark-Gluon Plasma (QGP) from hadronic matter, 
namely the deconfinement of color charges and the restoration of the spon- 
taneously broken chiral symmetry (which is believed to generate most of 
the visible mass in the universe). 

Heavy-quark (HQ) observables are expected to provide new insights 
into the aforementioned problems (with the exception of chiral symme- 
try restoration). On the one hand, charm (c) and bottom (6) quarks, due 
to their relatively large mass, should suffer less energy loss and undergo 
delayed thermalization when traversing the QGP, and thus be more sensi- 
tive to the interactions with thermal partons than light quarks and gluons. 
Therefore, it came as a surprise when the measurement of single-electron 
(e^) spectra associated with the semileptonic decay of open-charm (and 
-bottom) hadrons showed a factor 4-5 suppression 4,5 , very comparable to 
what has been found for pions. In addition, the observed values for the el- 
liptic flow, uf (Pt)j reach up to 10% at pr— 2 GeV, indicating the build-up of 
substantial early collectivity of c quarks. These observations reinforced the 
implementation of elastic energy-loss processes into the theoretical descrip- 
tion of the spectra 6-10 . On the other hand, HQ bound states (quarkonia), 
due to their small size/large binding, are suitable probes of a surrounding 
medium of sufficiently high density. E.g., for a typical (ground-state) char- 
monium of size r=0.25fm, the relevant parton density at which significant 
modifications are to be expected is n«3/(47rr 3 )« 10-20 fm -3 , which for an 
ideal QGP (with Nf =2.5 massless flavors) translates into a temperature of 
T«270MeVwl.5T c . Indeed, recent lattice QCD (1QCD) calculations indi- 
cate that J/ip and r\ c states survive in a QGP up to ~2T C . The theoretical 
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challenge is then to disentangle and quantify (a) medium modifications of 
the binding potential (e.g., color screening), (b) parton-induccd dissocia- 
tion reactions, and (c) in-medium changes of the HQ mass which affects 
both the bound-state mass and its decay threshold, as well as to identify 
and establish connections to observables in heavy-ion collisions. 

In the first part of the talk (Sec. 2), we evaluate HQ diffusion and ther- 
malization in the QGP employing both pcrturbative and nonperturbative 
elastic interactions. This problem is particularly suited to a Fokker-Planck 
approach for Brownian motion in a thermal background 11 , which allows 
to describe both the quasi-thermal and kinetic regime, and its transition. 
A Langevin simulation for c and b quarks in an expanding QGP at RHIC 
is supplemented with coalescence (and fragmentation) at T c , thus imple- 
menting all of the three main features (i)-(iii) of the initial RHIC data. 
In the second part of the talk (Sec. 3), we discuss in-medium properties 
of heavy quarkonia, including information from 1QCD. Wc outline a T- 
matrix approach in which HQ potentials from 1QCD are used as input to 
simultaneously describe bound and scattering states. This enables a com- 
prehensive evaluation of euclidean correlation functions, which in turn can 
be checked against rather accurate results from 1QCD. We briefly discuss 
how in-medium quarkonium properties reflect themselves in observables in 
ultrarclativistic heavy-ion collisions. Sec. 4 contains our conclusions. 

2. Open Charm and Bottom in the QGP 

Recent calculations of hadronic spectral functions in a QGP, both within 
1QCD and lQCD-based potential models, suggest that resonance (or bound) 
states persist up to 2T C , for both heavy- (QQ) and light-quark (qq) systems. 
We conjectured 6 that this also holds for heavy-light systems, and therefore 
could lead to significantly faster thcrmalization for c and b quarks as com- 
pared to perturbative QCD (pQCD) processes. Even at T=350MeV (which 
roughly corresponds to initial temperatures at RHIC) and for a s =0A, clas- 
tic pQCD scattering, which is dominated by i-channel gluon exchange, re- 
sults in a thermal relaxation time for c quarks above lOfm/c, well above 
a the typical QGP lifetime of TQGP~5fm/c. Recent data on single-electron 
spectra associated with semileptonic heavy-meson decays have corroborated 
the need for nonperturbative HQ interactions in the QGP. As shown in 
Fig. 1, perturbative energy-loss calculations appear insufficient to describe 
the strong suppression in the pr-spectra (left panel, using both elastic and 
radiative interactions) and the large V2{pt) (right panel, using radiative 
energy loss with a substantially upscaled transport coefficient). While the 
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Fig. 1 
heavy- 



Pcrturbativc QCD energy-loss calculations for spectra from semilcptonic 
meson decays. Left panel: nuclear suppression factor 10 , right panel: elliptic flow 12 . 



plots indicate a (lower) applicability limit of the pQCD energy-loss ap- 
proach of p|n~2 GeV, we will argue below that nonperturbative effects may 
be relevant (or even dominant) up to at least pfn~5GeV. In particular, 
radiative energy-loss calculations typically do not include the backward 
reactions (detailed balance) which are essential for building up collective 
behavior of the heavy quarks in the expanding thermal medium. As men- 
tioned above, in this context a Brownian-motion approach is well suited 
to address the thermalization of the heavy quarks 6,7 ' 11 ' 13 . Upon expand- 
ing the Boltzmann equation in small momentum transfers, one arrives at a 
Fokker-Planck equation for the HQ distribution function, /, 

d£ _ d(pf) (P£ 

dt" 1 dp + p cV ' ( j 

with (momentum) drag (7) and diffusion (D p ) constants (related via T = 
D p /jMq). The latter are calculated from corresponding matrix elements 
for HQ scattering off light partons. As mentioned above, our main ingredient 
here are resonance-mediated elastic interactions in s- and it-channel (cf. left 
panel of Fig. 2) which we model with an effective Lagrangian according to 6 

C = -G QT ^-^ $ q + h.c ; (2) 

$ is the resonance field representing and 11 B" mesons, which we dress 
by evaluating the corresponding one-loop (Q-q) selfenergy. Parameters of 
the model are the resonance mass (which we fix at 0.5 GeV above the Q-q 
threshold) and coupling constant, G, which determines the resonance width. 
HQ and chiral symmetry (above T c ) imply degeneracy of the pseudo- / scalar 
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Fig. 2. Elastic HQ-scattering processes via D-(B-)meson like resonances (left panel) 
and resulting thermalization times compared to pQCD results (right panel) 6 . 

and axial-/vector channels represented by the Dirac matrices, T. The tem- 
perature dependence of the resulting thermal relaxation times, tq = 7 _1 , 
for c- and fo-quarks shows a factor of ^3 reduction compared to calculations 
where only pQCD elastic scattering is accounted for (right panel of Fig. 2). 
This is very significant for c-quarks as r c is now comparable to (or even 
below) the duration of the QGP phase at RHIC (tqgp~5 fm/c), whereas 
tji>tqgp still, due to the large fe-quark mass. 

The T- and p-dependent drag and diffusion coefficients have been imple- 
mented into a relativistic Langevin simulation 6 using an (elliptic) thermal 
fireball expansion for (semi-) central Au-Au collisions at RHIC as a back- 
ground medium, where the bulk flow is adjusted to hydrodynamic simula- 
tions in accordance with experiment. At the end of the QGP phase (which 
terminates in a mixed phase at T c =180MeV), the c- and &-quark output dis- 
tributions from the Langevin simulation are subjected to hadronization into 
D- and B-mesons using the quark-coalescence model of Ref. 14 . Since the 
probability for coalescence is proportional to the light-quark distribution 
functions, it preferentially occurs at lower pt; "left-over" heavy quarks are 
hadronized via <5-function fragmentation. The extra contribution to momen- 
tum and i>2 from the light quarks increases both the Raa{pt) and V2(pt) 
of the heavy-meson spectra, relative to a scheme with fragmentation only. 
Rescattering in the hadronic phase has been neglected. The resulting D- and 
i?-meson spectra are decayed semileptonically resulting in a nuclear sup- 
pression factor and elliptic flow which compare reasonably well with recent 
RHIC data 4,5,15 up to pr— 5 GeV (Fig. 3). At higher momenta, radiative 
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Fig. 3. Calculations for single-electron spectra arising from scmileptonic heavy-meson 
decays. Left panel: Nuclear suppression factor right panel: elliptic flow. 



energy loss (not included here) is expected to contribute significantly 51 . 

Another issue in interpreting the scmileptonic e ± spectra concerns the 
relative contributions of charm and bottom. Since bottom quarks are much 
less affected by rescattering effects, their contribution reduces the effects 
in both Raa and «2 appreciably. Following Rcfs. 9 ' 17 , we have deter- 
mined the initial charm spectra by reproducing D- and D*-spectra in d- 
Au collisions 18 , calculated the pertinent e ± spectra and assigned the miss- 
ing yield relative to the data 19 to bottom contributions. In this way, the 
crossing between charm and bottom contribution in p-p collisions occurs 
at pt—5 GeV. Due to the strong quenching of the c-quark spectra in the 
QGP this crossing is shifted down to pr— 2- 5-3 GeV in Au-Au collisions. 
For charm contributions only, our nonperturbative rescattering mechanism 
results in a pertinent electron Raa as low as 0.1 and an elliptic flow as 
large as 14%. Thus, the severeness of the u e ± puzzle" partially hinges on 
the baseline spectra in more elementary systems (p-p and d-Au). Obviously, 
a direct measurement of D-mesons in Au-Au would be very valuable. 



3. Heavy Quarkonia in the QGP 

3.1. Potential Models, Spectral Functions and Lattice QCD 

Heavy quarkonia (c-c and b-b bound states) have long been identified as suit- 
able objects to quantitatively investigate the properties of QCD, including 
color confinement 20 . A particularly appealing feature is the applicability 
of potential models as a nonrelativistic effective-theory approximation to 



a Note that in a chemically equilibrated QGP, the effects of induced gluon radiation are 
presumably smaller than in a pure gluon plasma as assumed in Refs. 10 > 12 ; an estimate 
of (some) 3-body scattering diagrams has been performed in Ref. 16 . 
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QCD (cf. Ref. 21 for a recent review). In the vacuum, the heavy-quarkonium 
spectrum is indeed reasonably well described by the "Cornell potential" , 
consisting of a Coulomb plus linear confining term. Subsequently it was re- 
alized that, when immersing quarkonia into a medium of deconfincd color 
charges, Debye screening will reduce the binding and eventually lead to the 
dissolution of the bound state 22 ' 23 ; already slightly above T c , a substantial 
(factor ~2-3) reduction in charmonium and bottomonium binding energies, 
£b, has been found, together with the possibility that the (5-wave) ground 
states (r/ c ,6, J/ip, T(15)) survive (well) above T c . Modern 1QCD (-based) 
calculations qualitatively support this picture: On the one hand, directly 
extracted charmonium spectral functions (in quenched approximation) in- 
dicate resonance peaks up to ~2T C 24 ' 25 . On the other hand, (both quenched 
and unquenched) 1QCD results for the (color-singlet) free energy, 

F QQ (r;T) = U QQ (r;T)-TS QQ (r;T) , (3) 

have been implemented into potential models within a Schrodinger 26-29 or 
Lippmann-Schwinger equation 30 . There is an ongoing debate as to which 
quantity {Fqq, the internal energy Uqq, or a linear combination thereof 31 ) 
is the most appropriate one to use as a potential. Maximal binding is ob- 
tained with Uqq, and even in this case, the ground-state charmonium (bot- 
tomonium) binding energies are reduced to Eb — 0.1-0. 2 GeV (0.3-0.6 GeV) 
at T ~1.5 T c , compared to ~0.6GcV (1.1 GeV) in vacuum. However, a 
quantitative description of in-medium modifications of quarkonia requires 
at least two additional components, i.e., the effects of parton-induced dis- 
sociation reactions and of in-medium HQ masses. The objective is thus to 
construct an approach that (i) is consistent with finite-T 1QCD and (ii) 
allows for reliable applications to observables in heavy-ion collisions. Only 
then the original idea 22 of using quarkonia as a tool to identify (and char- 
acterize) confinement (or at least color screening) may be realized. 

The interplay of medium effects becomes more transparent in terms of 
the charmonium propagator, schematically written as (nonrelativistic) b 

D*(w;T) = [w-(2f7£- en) +iT*] _1 . (4) 

It shows that screening (affecting the binding energy and thus the real part 
of the inverse propagator), and parton-induced dissociation (governing the 
width, r^, of the spectral function, Im D^) arc not mutually exclusive but 
have to be taken into account simultaneously 32 . 



b Writing Eq. (4) in this form implies that the Q-Q potential is strong enough to generate 
a pole in the scattering amplitude. A microscopic (T-matrix) approach is discussed below. 
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Fig. 4. Lifetimes of J/ip (left panel, gluo-dissociation vs. quasifree dissociation 
mechanisms 34 ) and Y states (right panel, quasifree dissociation 35 ) in the QGP. 



The inelastic width (or dissociation rate) can be expressed as 

r* = (r*)- 1 = J ^ /*( Wfcj T) v rci 4 iss (s) (5) 

in terms of (thermal) parton-distribution functions, f p (p~q,q, g), and the 
parton-induced break-up cross section, a^, lss . The latter has first been eval- 
uated for gluo-dissociation 33 , g + *f? — » c + c, using Coulomb wave- functions 
for the quarkonia, leading to a T-dependent J/ip lifetime as shown in the 
left panel of Fig. 4. For the vacuum J/ip binding energy, £#=0.63 GeV, the 
lifetime is well behaved showing a rather steep decrease with T (dashed 
line). However, if an in- medium binding energy is used as following from 
screened (or 1QCD) potentials, the lifetime is initially reduced (dotted line), 
but if the binding energy drops below ~0.1 GeV, the phase space in the gluo- 
dissociation reaction strongly shrinks leading to an unphysical increase of 
the J/ip lifetime with T. This artifact is even more pronounced for the less 
bound excited states (% c , ?/;'). This problem has been remedied by intro- 
ducing the "quasifree" dissociation process 34 , p + ^—tp + c+c, where a 
partem knocks out the Q (or Q) from the bound state. With reasonable val- 
ues for the strong coupling constant (a s ~0.25) the resulting charmonium 
widths are in the ~0.1GeV range for T~1.5 T c (as relevant for RHIC) 
with a well-behaved T-dependence. The order of magnitude of the inelastic 
width is easily reproduced using a simplified estimate of Eq. (5) according 
to r^wn p (T)adi ss i'rel) with a parton density n p ~10fm -3 (at T«0.25GeV), 
a cross section of 1 mb and u re i=l/2 one finds r*«0.1 GeV. 

When applied to bottomonia, the same arguments apply 35 ; the right 
panel of Fig. 4 illustrates the impact of in-medium binding energies on bot- 
tomonium dissociation: in the RHIC temperature regime, screening may 
lead to up to a factor of ~20 (!) reduction in the ground-state T(15) life- 
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Fig. 5. Finite-temperature euclidean correlation functions computed in 1QCD in the 
pseudoscalar (rj c , left panel) and scalar (xc right panel) channels 25 . 

time, bringing it down to ^2-5fm/c at T=0.3-0.4GeV. This has impor- 
tant consequences for bottomonium suppression at RHIC (as elaborated in 
Sec. 3.2 below), rendering it a sensitive probe of color screening 35 . 

A key issue for model descriptions of quarkonium spectral functions in 
the QGP is the implementation of constraints from 1QCD. The standard 
object computed on the lattice is (the thermal expectation value of) a 
hadronic current-current correlation function in euclidean time, r (and at 
zero 3- momentum), in quantum- number channel, a 36 , 

G a (r;T) = (j a (T)jl(0)) . (6) 

The connection to the spectral function, a a , is given by a convolution with 
a temperature kernel as 

cosh[w(r - 1/2T)] 



G a (t; T) 



da; er Q (w; T) 



(7) 



sinh[w/2T] 

While euclidean correlators can be computed with good accuracy (Fig. 5), 
their limited temporal extent at finite T, 0<r<l/T, renders the extraction 
of (Minkowski) spectral functions more difficult. On the contrary, model 
calculations of spectral functions are readily integrated via Eq. (7) and 
compared to the rather precise euclidean correlators from 1QCD 37 - 29 . To 
better exhibit the T-dependence induced by the spectral function, both 
lattice and model results for G a (r;T) are commonly normalized to the so- 
called "reconstructed" correlator, which follows from Eq. (7) by inserting a 
free spectral function, a a (uj,T = 0). 

Eq. (6) requires the knowledge of the quarkonium spectral function at 
all (positive) energies. A schematic decomposition consists of a bound-state 
part and a continuum, 



2M l Ff5{uj 2 - M?) + ^ lu 2 / cont (w, E thY )e(uj - E thl ) , (8) 
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where, for fixed a, i runs over all bound states with masses Mj and couplings 
Fi, and the threshold energy (Pthr) an d function / CO nt characterize the onset 
and plateau of the continuum (usually taken from pQCD). 

In Rcf. 29 , T-dependent bound-state masses and couplings have been 
evaluated by solving a Schrodingcr equation with in-medium potentials 
(screened Cornell or 1QCD internal energy, Uqq); the threshold energy 
has been inferred from the asymptotic (r ~ > oo) value of the respective 
potential, E^ d (T) = 2m c + V(r — > oo; T) (which monotonously decreases 
for T > T c ). For the scalar (% c ) channel it has been found that, despite a 
rather quickly dissolving bound-state contribution, the correlator increases 
significantly (consistent with the lattice result in Fig. 5, right panel), due to 
the decreasing continuum threshold. The latter also leads to an increasing 
pscudoscalar correlator, which is not favored by 1QCD (Fig. 5, left panel). 

In Ref. 38 the separation of bound-state and continuum parts is im- 
proved by employing a T-matrix approach 39 ' 30 for the Q-Q interaction, 



T a (E;q',q) = V a (q',q)-- / dk k 2 V a (q' , k) G QQ (E; k) T a (E; k, q) . (9) 



For V a the Fourier transformed and partial-wave expanded 1QCD internal 
energy has been used. Spin-spin (hyperfine) interactions are neglected im- 
plying degeneracy of states with fixed angular momentum (S'-wave: rj c and 
J/tp, P-wave: Xc0,i,2)- Gqq(E; k) denotes the intermediate 2-particle prop- 
agator including quark selfencrgics. The correlation and spectral functions 
follow from closing the external legs with 3-momcntum integrations as 



G a (E) = / Gqq + / Gqq / T a G QQ , a a (E) = a lmG a {E) , (10) 



where the coefficient a depends on the channel a. For a fixed c-quark mass 
of m c =1.7GeV, the resulting charmonium spectral functions confirm that 
bound states are supported in the S'-wave up to ~3T C while dissolved in the 
P-wave below 1.5 T c , see Fig. 6. In both cases, however, one finds a large 
(nonperturbative) enhancement of strength in the threshold region. The 
corresponding euclidean correlators are displayed in Fig. 7, normalized to a 
reconstructed one using the vacuum bound-state spectrum and a perturba- 
tive continuum with onset at the free DD threshold, E^=2mD— 3. 75 GeV. 
The Xc correlator (right panel) is well above one, both due to a lower in- 
medium threshold, P^ d =2m c =3.4 GeV, and the nonperturbative rescat- 
tering enhancement. While this is qualitatively consistent with 1QCD, the 
temperature dependence is opposite, indicative for a further threshold de- 
crease with temperature. In the r\ c channel (left panel), the significant re- 
duction in binding energy (due to color screening) reduces the correlator 




CO 
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Fig. 6. S- (left panel) and P-wave (right panel) charmonium-spectral functions from 
the T-matrix approach with an lQCD-internal energy as potential 38 . 
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Fig. 7. Normalized cuclidcan correlation function for S- (left panel) and P-wave (right 
panel) charmonium states within a T-matrix approach using lQCD-internal energies 38 . 



at large r, qualitatively in line with 1QCD. In fact, the latter shows less 
reduction, leaving room for a T-dependent reduction of E™ d as well. 

The scattering equation (9) is furthermore well suited to study finite- 
width effects. E.g., when dressing the c-quarks with an imaginary self- 
energy corresponding to r c =50 MeV, inducing a charmonium width of 
r^~100MeV (as in the left panel of Fig. 4), the euclidean correlators vary 
by no more than 2%. This indicates that there is rather little sensitivity to 
phcnomcnologically relevant values for the quarkonium widths. 



3.2. Quarkonium Phenomenology in Heavy-Ion Collisions 

The formation of thermalized matter above T c in ultrarelativistic heavy- 
ion collisions (recall the discussion in the introduction) provides the basis 
for describing the production systematics of quarkonia in terms of their 
in-medium properties as discussed in the previous section. The suppression 
of the initially produced number of quarkonia, N^,, may be schematically 
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Fig. 8. Centrality dependence of the nuclear modification factor for J/ip (left panel) 40,46 
and T (right panel) 35 at RHIC. 



written as proceeding through 3 stages, 

S to t = exp(-a^l g N L) exp(-r QG p t qgp ) exp(-r H c Thg)- (H) 

With the suppression in the hadron gas (HG) believed to be small, and the 
( "pre-cquilibrium" ) nuclear absorption (cross section) inferred from p- A 
experiments, the suppression factor directly probes the dissociation rate of 
^'s in the QGP. The observed J/^-suppression pattern at SPS can indeed 
be well described using the quasifree dissociation process with in-medium 
binding energies consistent with lQCD-potcntial models 34,40 . 

At RHIC energy, copious production of cc pairs (iV c g~20 in central 
v / i=200 AGeV Au-Au 41 , compared to 0.2 at SPS) opens the possibility 
of secondary charmonium production 42,34 . This is, after all, required by 
detailed balance in the dissociation reaction, "f+.g ^ c+c+X, provided 
the \& state still exists at the given temperature 40,43 . If, in addition, the 
c-quarks are close to thermal equilibrium, one can apply the rate equation, 

dN* 



dt 



K q (T)] 



(12) 



for the time evolution of the J/ip number. Pertinent predictions (upper solid 
line in the left panel of Fig. 8) agree reasonably well with current PHENIX 
data 44 suggesting that the regeneration component becomes substantial in 
central Au-Au. Corrections due to incomplete c-quark thermalization, as 
well as a lower dissociation (and thus formation) temperature, reduce the 
regeneration (lower solid curves in the left panel of Fig. 8). It has also 
been suggested 45 that available J/tp data at SPS and RHIC are compatible 
with suppression of only Xc and ip' (which make up ~40% of the inclusive 
J/ip yield in p-p), with the direct J/-0's being unaffected. However, if at 
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RHIC the J/tp width is of order 0.1 GcV (cf. previous section) and the QGP 
lifetime at least 2fm/c, one finds a QGP suppression factor of Sqgp— 0.4. 

Finally, turning to T production, the increase of the T dissociation rate 
due to color screening (recall right panel of Fig. 4) could lead to a T yield 
in central Au-Au (right panel of Fig. 8) which is more suppressed for the 
J ftp (left panel), and thus provide a rather striking QGP signature 35 . 

4. Conclusions 

Charm and bottom hadrons are valuable sensors of the sQGP. In the open- 
flavor sector, the inadequacy of pQCD to account for the suppression and 
collectivity of current e spectra at RHIC may be a rather direct indica- 
tion of nonperturbative rescattering processes above T c . We have elaborated 
these in a concrete example using (elastic) resonance interactions. The lat- 
ter will have to be (i) evaluated more microscopically, e.g., using input 
interactions from 1QCD, and (ii) combined with radiative energy loss. 

In the quarkonium sector, steady progress is made in implementing first- 
principle information from finite-T lattice QCD into effective models suit- 
able for tests in heavy-ion reactions. T-matrix approaches incorporate ef- 
fects of color-screening, parton-induced dissociation, and in-medium masses 
and widths of heavy quarks (which, in turn, connect to the open-flavor sec- 
tor). We have suggested an sQGP signature in terms of stronger suppression 
for T relative to J/ip, which would be a direct proof of J/ip regeneration 
and carries large sensitivity to color screening of bottomonia. 
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